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II.  WEDGES  AND  CONES 

by 

* 

Michael  L.  Shorenstein 

Abstract 

An  analysis  is  given  for  the  viscous  hypersonic  flow  past  slender 
sharp  wedges  and  unyawed  cones.  The  loading  edge  or  nose  region 
treated  encompasses  the  merged  layer  regime  in  which  the  shock  thick¬ 
ness  is  not  small  in  comparison  with  the  viscous  layer  thickness  but 
for  which  in  the  axisymmetric  case  the  shock  thickness  is  small  in 
comparison  with  the  radius  of  transverse  shock  curvature.  A  Navier- 
Stokes  flow  is  assumed,  and  both  the  shock  structure  and  viscous 
layer  are  taken  to  be  locally-similar  of  boundary  layer  type.  The 
analysis  follows  that  used  by  Shorenstein  and  Probstein  for  the 
corresponding  flat  plate  problem.  For  axisymmetric  cone  flow,  a 
modification  of  a  transformation  introduced  by  Probstein  and  Elliott 
is  given  to  reduce  the  axisymmetric  viscous  layer  equations  includ¬ 
ing  transverse  curvature  to  ’exactly"  twc  -dimensional  form.  Results 
for  the  state  conditions  behind  the  shock,  and  the  wall  pressure, 
skin  friction,  and  heat  transfer  rate  are  found  to  be  in  good  agree¬ 
ment  with  available  experimental  data. 

*  The  author  wishes  to  thank  sincerely  Professor  Ronald  F.  Probstein 
of  M.I.T.  for  his  support  during  the  course  of  this  investigation 
and  for  his  valuable  suggestions  and  criticisms.  This  research  was 
supported  by  the  Advanced  Research  Projects  Agency  and  was  tech¬ 
nically  administered  by  the  Office  of  Naval  Research  under  Contract 
No.  N0001 4-0204-0040,  ARPA  Order  No.  322. 

Research  Assistant,  M.I.T. ,  presently  at  Northern  Research  &  Engin¬ 
eering  Corp.,  Cambridge,  Mass. 
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Nomenclature 

■  velocity  ratio,  u*/U 

8  °® 

■  given  by  Eq.  (18) 

-  Chapman-Rub es in  constant 

■  given  by  Eq.  (20a) 

■  given  by  Eq.  (20b) 

■  skin  friction  coefficient 

-  specific  heat  at  constant  pressure 

■  reduced  density,  p/p 

00 

-  8 treamf unction  given  by  Eq.  (6) 

■  reduced  slip  velocity,  u*  /u* 

sfc  s 

-  given  by  Eq.  (43) 

-  shear  stress,  Nf 

nn 

■  value  of  g  at  body  surface 
«  specific  enthalpy,  cpT 

-  total  enthalpy 

-  given  by  Eq.  (11) 

■  index  taken  to  be  zero  for  two-dimensional  flow  and  unity 

for  axl8ymmetric  flow 

■  arbitrary  fixed  reference  length 

■  frees tr earn  Mach  number 

■  normal  coordinate  for  shock  structure 

■  given  by  Eq.  (5) 

■  static  pressure 
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■  reduced  pressure, 

■  Prandtl  number 

-  given  by  Eq.  (28) 

■  given  by  Eq.  (29) 

*  given  by  Eq.  (42) 

■  radius  of  transverse  curvature,  r,  +  y*  cos  9, 

d  b 

■  gas  constant 

■  local  Reynolds  number,  p  U  x*/p 

CD  00  00 

-  Stanton  number 

-  static  temperature 

-  freestream  stagnation  temperature 

-  velocity  component  parallel  to  shock  surface  and  to  body 

surface,  respectively  [see  Fig.  2] 

■  freestream  velocity 

■  velocity  component  normal  to  shock  surface  and  to  body 

surface,  respectively  [see  Fig.  2] 

-  transformed  normal  velocity  given  by  Eq.  (2a) 

■  reduced  velocity,  v/U^ 

■  reduced  velocity,  u/Uw 

■  coordinate  along  body  surface  and  normal  to  body  surface, 

respectively  [see  Fig.  2] 

■  transformed  coordinates  given  by  Eq.  (2a) 

■  Inverse  Reynolds  number  (p  U  r,  /p  )  ^ 

00  00  Q  00 

-  specific  heat  ratio 

■  shock  thickness 

■  stretched  normal  coordinate  for  shock  structure,  Eq.  (23) 

■  stretched  normal  coordinate  for  viscous  layer,  Eq.  (6) 
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transformed  location  of  the  body  surface  for  the  case  of  slip 
location  of  the  shock  surface 

body  surface  Inclination  angle  with  respect  to  the  free 
stream 

shock  surface  inclination  angle  with  respect  to  the  body 
surface 

reduced  enthalpy  given  by  Eq.  (6) 
viscosity 

given  by  transformation  Eq.  (^b) 

transformed  coordinate  along  body  surface  given  by  Eq.  (6) 
density 

cross-section  of  streamtube  incident  to  viscous  layer 

reduced  enthalpy  given  by  Eq.  (23) 

3  1/2 

interaction  parameter,  M  (C/Re  .  ) 

00  x* ,  “ 


Subscripts 


b 

nsfc 

s 

si 

n 


condition  at  the  body  surface 
quantity  obtained  from  no-slip  solution 
condition  at  the  shock  surface 
quantity  obtained  from  slip  solution 
derivative  with  respect  to  n 
condition  at  the  freestream 
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I.  Introduction 

Two  classes  of  continuum  flow  models  have  been  investigated 
recently  for  treating  the  hypersonic  flow  past  sharp  bodies  up¬ 
stream  of  the  strong  Interaction  regime  In  the  merged  layer  [see 

Fig.  1]  where  the  shock  thickness  must  be  considered  finite  In 

1-3 

comparison  with  the  viscous  layer  thickness.  One  class  applies 

the  full  Navler-Stokes  equations  In  finite  difference  form  In  order 

to  determine  to  what  extent  a  continuum  may  meaningfully  describe 

the  flow  field  In  the  merged  layer.  This  complex  scheme  requires 

rather  large  computer  storage  space. 

Another  class  of  merged  layer  models  uses  the  Navler-Stokes 

equations  but  assumes  "local  similarity"  to  apply  in  both  the 

viscous  layer  and  shock  wave,  that  Is,  streamwlse  variations  of  all 

properties  are  taken  to  be  sufficiently  small,  so  that  the  Navler- 

Stokes  equations,  when  applied  locally,  can  be  approximated  by 

ordinary  differential  equations.  This  description  of  the  merged 

viscous  shock  layer  Is  obtained  from  locally-slmllar  solutions 

"patched"  at  the  Interface  between  the  outer  shock  structure  portion 

and  the  Inner  viscous  layer  portion  in  which  the  equations  are  taken 

to  be  of  boundary  layer  type  with  normal  pressure  gradient  effects 

neglected.  Results  In  Ref.  6  from  such  a  merged  layer  model  for  a 

7-9 

flat  plate  have  compared  well  with  experimental  data  .  The 
theoretical  calculations  have  shown  cne  relative  Importance  of 
several  rarefaction  effects  over  a  cold  sharp  flat  plate  aligned 


parallel  to  the  hypersonic  freestream.  Dominant  among  these  effects 
was  shown  to  be  the  finite  shear  at  the  Interface  between  the  back 
of  the  shock  and  the  outer  edge  of  the  viscous  layer.  The  effect  of 
Including  velocity  slip  and  temperature  jump  at  the  body  surface  was 
shown  to  be  smaller  but  non-negligible,  while  the  effect  of  including 
longitudinal  shock  curvature  was  found  in  comparison  to  be  negligible. 

The  present  analysis  extends  this  flat  plate  continuum  merged 
layer  model  of  Ref.  6  to  wedge  and  axisymmetrlc  conical  flows  by  in¬ 
vestigating  the  additional  effects  of  body  slope  and,  for  conical 
flow,  the  effects  of  axial  symmetry  and  transverse  curvature  effects. 
Longitudinal  shock  curvature  effects  are  neglected  as  small  based  on 
the  results  of  Ref.  6.  Body  slope  effects  are  characterized  by  the 
parameter  6b<  where  6^  is  the  half  angle  of  the  wedge  or  unyawed 
cone.  The  analysis  is  extended  only  to  slender  bodies  with  9fe  <<  1 
for  which  the  normal  component  of  the  viscous  layer  velocity  can  be 

neglected  compared  to  the  tangential  component.  This  restriction  is 

• 

introduced  in  order  to  solve  the  viscous  layer  problem  by  the  method 
of  Ref.  6. 

In  the  axial  symmetric  case,  a  modification  of  a  dangler  type 
transformation  given  by  Probstein  and  Elliott10  is  introduced  to 
reduce  the  axisymmetrlc  viscous  layer  equations  to  an  "exactly"  two- 
dimensional  form.  The  modification  Involves  a  transformation  applied 
to  the  viscosity  coefficient,  and  this  transformation  is  shown  to 
include  the  effects  of  transverse  curvature  characterized  by  the 
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parameter  M^a.  Here  a  is  an  inverse  Reynolds  number  based  on  the  body 
radius  r,  of  transverse  curvature,  and  M  is  the  freestream  Mach 

D  oo 

number.  A  ijew  term  proportional  to  M  o  is  found  to  modify  the  reduced 

00 

viscous  layer  momentum  equation  of  Ref.  6  and  is  shown  to  account  for 
transverse  curvature  effects. 

For  the  present  investigation,  two  regimes  of  M  o  are  dis- 

00 

tinguished:  I,  a  regime  immediately  adjacent  to  the  sharp  nose  where 
M  a  >  1  and  II,  a  regime  downstream  of  the  sharp  nose  where  M  a  <<  1. 

00  »  oo 

Only  the  second  regime  is  treated  here.  This  restriction  is  introduced 

to  obtain  a  simple  closed  form  shock  structure  solution  by  a  method 

analogous  to  Ref.  6.  In  the  conical  shock  structure  equations, 

transverse  shock  curvature  terms  are  of  order  A  / r  where  A  is  the 

s  s  s 

shock  thickness  and  r  is  the  radius  of  transverse  shock  curvature. 

s 

It  is  found  that  if  M  a  «  1,  then  terms  of  order  A  /r  are 

"  s  s 

negligible  compared  to  the  other  terms  in  the  shock  structure 
equations.  These  equations  can  then  be  reduced  to  the  form  of  the 
"zero-order"  t:hock  equations  solved  in  Ref.  6. 

II.  Viscous  Layer 

As  in  Ref.  6,  the  merged  shock  layer  1?  described  by  a  viscous 
layer  of  boundary  layer  type  patched  at  an  outer  interface  (the  "shock 
surface")  to  a  thick  shock  structure.  The  velocity  and  state 
variables  from  the  solution  of  the  viscous  layer  equations  are  matched 
locally  at  the  shock  surface  to  the  corresponding  quantities  from  the 
solution  of  the  shock  structure  equations. 


The  governing  equations  for  the  viscous  layer  written  in  the  local 
body-centered  coordinate  frame  (x*,y*)  depicted  in  Fig.  2  are 


fe*  <rJeu*>  +  I7*  (rjpv*)  -  0 


,(u*  +  v*  _  d£_  9_  ,  3c* .  ji  3t_  3u* 

,vu  av*  v  a«*'  dx*  gy*  gy*'  J  »•  a«*  a»* 


r  3y*  3y* 


p(u*  Ql.  +  v*  1L_\  .  u*  d£_  +  L_  rl —  /u 

3x*  3y*  dx*  Pr  3y*  'y  3y*' 


+  1  1L  .52L_  +. 

J  r  3y*  3y*J  y'‘3y*; 


where  j  ■  0  for  two-dimensional  flow  and  j  *  1  for  axisymmetric  flow, 
with  r(x*,y*)  ■  rfe(x*)  +  y*  cos  efa  the  local  radius  of  transverse 
curvature  for  the  case  of  axisymmetric  flow  with  0b  the  body  slope 
relative  to  the  freestream  direction.  The  tangential  and  normal  com¬ 
ponents  of  velocity  are  u*  and  v*,  respectively  [see  Fig.  2],  h  is 
the  specific  enthalpy,  p  is  the  pressure,  pis  the  viscosity,  and  Pr 
is  the  Prandtl  number  taken  to  be  a  constant. 

Probstein  and  Elliott^  introduced  a  generalized  Mangier  trans¬ 
formation  which  reduced  Eqs.  (1)  to  an  almost  two-dimensional  form 
by  applying  the  transformations 


dx  -  (j—)2^  Jx*  , 


3x  3y 


dy  ■  (j^  dy* 


L  3y 


V  -  (~ 2>j  V*  +  j((^-)2  |jyl  u* 

r,  b 
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where  L  is  an  arbitrary  fixed  reference  length.  By  also  introducing 
the  transformation 

_  2jL  y  cos  8. 

p  -  (1  +  - 2 - “)  U  (2b) 

rb 

Eqs.  (1)  can  be  reduced  "exactly"  to  two-dimensional  form.  This 

follows  from  the  fact  that  6.  and  r.  are  functions  of  x  only,  so 

b  b 

that  Eqs.  (1)  when  transformed  by  means  of  Eqs.  (2)  give  the 
homogeneous  system 

~  (pu*)  +  (pv)  -  0  (3a) 

3x  3y 


/  *  3u*  .  ~  3u*\  3  3u*  ,  dp  _ 

p(u*  +  v  — — )  -  —  (p  -— )  +  _  “  0  (3b) 

3x  3y  3y  3y  dx 

,  *  3h  —  3h.  13,—  3h.  —  ,3u*x2  .dp  _  . 

p(u*  —  +  v  — )  -  p-  —  (p  — )  -  p(~3j-)  -  u*  -  0  (3c) 

3x  3y  3y  3y  3y  dx 

which  are  "two-dimensional"  in  form.  It  must  be  pointed  out,  how¬ 
ever,  that  the  parameter  p,  although  having  the  dimensions  of  a 
viscosity,  is  not  a  true  state  variable  since  it  depends  not  only 
upon  temperature  but  also  on  the  geometry  of  the  system.  Still, 

Eqs.  (3)  can  be  interpreted  as  mathematically  analogous  to  the 
equations  of  two-dimensional  flow.  Furthermore,  Eqs.  (3)  can  be 
solved  by  the  method  used  in  the  viscous  layer  problem  of  Ref.  6 
if  v  can  be  taken  as  small  compared  to  u*.  The  present  work 


S-88K3WS9 
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therefore  extends  the  merged  layer  model  of  Ref.  6  only  to  slender 
bodies  (efa  «  1)  for  which  the  viscous  layer  streamline  curvature 
3y*/3x*  can  be  ignored  and  for  which  v*  is  negligibly  small  com¬ 
pared  to  u*.  For  such  bodies,  by  Eq.  (2a),  v  <<  u*. 

With  the  governing  equations  in  two-dimensional  form  and  with 
v  <<  u*,  the  analysis  can  more  directly  parallel  that  of  Ref.  6. 
Following  this  reference,  the  viscous  layer  adjacent  to  the  body  is 
assumed  locally  similar,  and  for  a  calorically  perfect  gas  with  a 
Prandtl  number  of  unity,  Eqs.  (3b)  and  (3c)  then  reduce  to 


(Nf  )  +  ff 

nn  n  nn 


-  0 


(4a) 


(»®)„  +  ■  0 


(4b) 


where 


N  -  Pu/PbWb 


(5) 


and  where  the  dimensionless  independent  and  dependent  variables  are 

e(*>  ■  I  w; d* 

1  w  0  0 

(6) 

f(n)  *  f  -j-  dn  (h)(h)  ■  — - ~ 

6  s  ^-“b 

Here  H  is  the  total  enthalpy  in  the  barred  system  (with  v  <<  u*)  and 


rnimmmmmmmmm.  1 
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the  subscripts  "s"  and  "b"  denote  quantities  at  the  shock  surface  and 
body  surface,  respectively. 

Following  Ref.  6,  we  first  consider  the  problem  neglecting 
velocity  slip  and  temperature  jump  at  the  body  surface.  With  the 
body  temperature  constant,  the  boundary  conditions  are  then 


f  -  f  -  0 
n 


■  0  at  n  ■  0 


fn  “  1  ®“  1  at  n  -  ng 


(7) 


Consider  now  the  energy  equation  (Ab)  which  is  satisfied  by 
the  Crocco  integral  -  f^.  From  the  definition  of  total  enthalpy 
(with  v  <<  u*) 


T_ 

T. 


(1  -  f  )  [1  +  A2  (-°)  f  ] 
n  Tb  n 


(8) 


where 


A(x)  ■  u*/U  (9) 

s  oo 

with  U  the  freestream  velocity.  For  the  hypersonic  conditions  con- 
sldered,  the  freestream  stagnation  temperature  Tq  is  given  by 
Tn  -  U2/2c  . 

0  »  p 

Consider  next  the  viscous  layer  momentum  equation  (Aa) .  An 
expression  for  N  Is  sought  so  that  Eq .  (Aa)  can  be  expressed  in 
terms  of  the  streamf unction  f(rO  and  its  derivatives  and  in  terms 
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of  Che  local  boundary  conditions  at  the  shock  surface  and  the  body.  As 

4 

shown  In  Pan  and  Probsteln  ,  the  normal  pressure  gradient  across  the 
viscous  layer  Is  negligible  when  viscous  layer  streamline  curvature 
Is  negligible.  Thus,  by  p  ■  pRT,  p.  /p  ■  T/T.  so  that  from  the  defini- 

D  D 

tlon  of  n [Eq.  (6)]  It  follows  from  Eq.  (8)  that  at  any  given  x 


y  -  l(201/2/pbu*l  Kn) 


(10) 


where 


n 

I(n)  -  |  (1  -  fp)  [1  +  A2(TQ/Tb)  f^J  dn 
0 


(11) 


Hence,  from  Eq.  (10)  and  from  the  definition  of  p  [Eq.  (2)] 


p/p  ■  1  +  BI(n) 


(12) 


where  B  Is  a  function  of  the  boundary  conditions  at  n  “  0  and 

n  -  ng.  Since  N  -  (p/pb)(p/pb)  -  (Tb/T)(p/pb)  ■  (Tb/T) (p/pb) (p/p)  and 

1/2 

If  a  viscosity-temperature  relation  of  the  form  p  «  T  Is  assumed, 
then  by  Eq.  (12) 

N  -  (Tb/T)1/2  [1  +  BI(n)]  (13) 

To  complete  this  expression  for  N,  we  now  obtain  B  In  terms 
of  the  local  boundary  conditions  at  the  shock  surface  and  the  body. 
This  Is  done  by  first  applying  local  conservation  of  mass  flux 


CTWWS 
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across  Che  shock  and  Into  the  viscous  layer  by  use  of 

7s 

P  U  o(J)  -  (2itL)J  f  pu*  dy  -  (2irL)j  (2?)1/2  f(n  )  (14a) 

«o  eo  J  S  s 

0 

where  o(j)  Is  the  cross-section  of  the  Incident  streamtube  given 
by 

cr(j  *  0)  -  [1  +  4  (cot  6  -  tan  0,  )  sin  20,  ]  y*/cos  0,  (14b) 

2S  D  D  S  D 

o(j  -  1)  -  irr2  (14c) 

s 

Now,  from  the  definitions  of  x  [Eq.  (2)]  and  of  £(x)  [Eq. 

(6)],  It  can  be  shown  that 

2C  -  2ptubu*  [(1  -  2j/3)  (rb/L)2J]  x*  (15) 

where  r,  -  x*  sin  0,  [see  Fig,  2]  has  been  used.  Next,  dividing 

D  D 

Eq.  (14a)  by  Eq.  (15)  for  the  case  j  ■  1,  Eq.  (10)  becomes 

y  -  {[2ii.f(n  )/p  U  ]  Or./r  )2  (2x*/3L)}  I(n)  (16) 

o  s  ao  ao  D  8 

Finally,  neglecting  longitudinal  shock  curvature,  we  set 

y*/x*  ■  tan  0  so  that 
s  s 

r  /x*  -  (r.  +  y*  cos  0,  )/x*  -  [1  +  (tan  0  /tan  0.  ))  sin  0, 

8  D  S  D  S  D  D 


(17) 
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and  by  comparing  Eq.  (12)  Co  Che  deflnlcion  (2b)  of  u  by  Che  use  of 
Eqs.  (16)  and  (17),  we  find 

B  ■  (8jy. /3p  Can  0,  [1  +  (Can  0  /Can  0,  )]^}  f(n  )  a  (18) 
b  <*>  b  8  b  s 

Here  a  Is  Che  Inverse  Reynolds  number  (p  U  r ,/\i  )  based  on  free- 

00  00  Q  00 

scream  condldons  and  Che  local  radius  r.  (x*)  of  Cransverse  body 

D 

curvature.  From  Eq.  (18),  B  vanishes  for  Cwo-dlmensional  flow 

(j  -  0). 

From  Che  energy  equaClon  In  Che  form  of  Eq.  (8),  from 
1/2 

U  ^  T  ,  and  from  Eqs.  (13)  and  (lo)  which  CogeCher  relaCe  Che 
parameCer  N  Co  Che  s Cr earn func Cion  f(n),  Co  lCs  derlvaclves,  and  Co 
Che  local  boundary  condldons  on  Che  viscous  layer,  Che  viscous 
layer  momenCum  equation  (Aa)  now  reduces  Co 


(1  +  C.I (n) }  f 
_ £ _ OIL 


{u-ym-c^n 


+  ff 
l/2'n  nn 


(19) 


where 


C1  “  A  <T0/Tb) 


(20a) 


C2  -  (8j  “)1/2  f(ng)/3Can  0b  [1  +  (Can  0s/can  0fa)]2}  M^a  (20b) 

0 


EquaClon  (19)  differs  from  Che  reduced  viscous  layer  momenCum  equaClon 
of  Ref.  6  only  by  Che  appearance  of  Che  new  Cerm  C2l(n)«  This  term 


*t,<r  -  'f  *  v  -ft-'  *  -  •  .(UMnuMiaMH 
■  ■  r ,'T|*3\Wv"Xnp5!W»l 
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introduces  the  effects  of  transverse  curvature  characterized  by  the 

parameter  M  a,  and  it  vanishes  for  the  case  of  two-dimensional  flow. 

00 

For  the  case  of  axial  symmetric  flow,  only  the  regime  M  o  <<  1  is 

00 

treated  in  this  work.  The  reason  for  this  will  be  discussed  in  con¬ 
junction  with  the  shock  structure  analysis. 

Now,  if  and  were  to  be  treated  as  parameters,  then  for 

each  value  of  (f  )  greater  than  the  value  which  gives  the  boundary 

nn  d 

layer  solution  (ng  ■+■«•),  a  viscous  layer  solution  would  be  obtained 

with  finite  shear  (proportional  to  Nf  )  at  the  shock  surface  n  *  n 

nn  s 

where  f  *  u*/u*  ■  1.  The  shock  surface  coordinates  x*,  y*  corres- 
n  s  s 

ponding  to  ng  are  determined  by  1)  employing  Eq.  (14a)  which  states 
that  mass  flux  into  the  viscous  layer  between  the  body  and  shock 
surface  is  conserved  and  2)  by  matching  the  velocity  u*  and  its 
gradient  3u*/Dy*  obtained  with  finite  shear  at  the  outer  edge  of  the 
viscous  layer  to  the  respective  quantities  obtained  at  the  inner  edge 
of  the  shock  from  the  shock  structure  solution. 

III.  Shock  Structure 

Since  the  flow  at  the  shock  surface  has  been  assumed  locally 

similar  within  the  viscous  layer  model,  all  shock  structure  variables 

at  the  inner  edge  of  the  shock  are  single-valued  functions  of  n  . 

s 

The  shock  transition  equations  in  the  coordinate  frame  centered 
at  the  shock  surface  are  illustrated  by  the  continuity  and  tangential 


momentum  equations 


VT- onMiwn) 


h  "J<">  ■ 0 
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(21a) 


1_  3_ 
rj  3n 


(21b) 


where  n  Is  the  coordinate  In  the  direction  measured  normal  from  the 
shock  surface  to  the  freestream,  and  the  velocity  components  u  and 
v  are  taken  normal  to  and  parallel  to  n,  respectively  [see  Fig.  2]. 
based  on  the  results  of  Ref.  6,  the  effects  of  including  longitudinal 
shock  curvature  terms  in  Eqs.  (21)  have  been  neglected  as  small. 

Integrating  the  continuity  equation  (21a)  yields 


r^pv  ■  constant 


-rJ 


p  v 
s  s  s 


r^p  v 

CD  00  00 


(22a) 


Substituting  Eq.  (22a)  into  Eq.  (21b)  and  multiplying  by  (r/r  )2^, 

s 

the  tangential  momentum  equation  (21b)  becomes 


,r  .j  3u  ,r  ,]  3  ,  ,r  . j  3u, 
Dsvs(-)  to  *  <r>  to  tu(-'  to1 


(22b) 


To  solve  these  shock  structure  equations,  the  following  reduced, 
dimensionless  variables  are  Introduced: 


f  p.U»  r 
—  (- 
I  u  r 


V  dn 


W  -  — 

u 


D 


0  IT 


V  -  2- 
U 


(23) 


4>  - 


h  +  -  v2 
U2 


% 


| 


/ 


/ 
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The  dimensionless  forms  of  Eqs.  (22)  by  (23)  are  then 


DV  «  [l+j(4  /r)]DV 
ss  J  s  s  »  « 


(24a) 


(24b) 


where  A  =  (r  -  r  )  denotes  shock  thickness.  Equations  (24)  differ 

S  00  9 

from  the  "zero-order"  shock  transition  equations  of  Ref.  6  only  by 
the  appearance  of  the  transverse  shock  curvature  parameter  A  /r  . 

8  S 

If  however  M  a  <<  1,  Chen  A  /r  <<  1  and  a  simple  closed  form  solu- 
*88 

tlon  to  Eqs.  (24)  can  be  readily  obtained  analogous  to  the  "zero- 
order"  shock  solution  In  Ref.  6.  In  the  present  model,  M^a  <<  1  Is 
satisfied  and  transverse  shock  curvature  terms  are  neglected  as 
small  In  the  shock  structure  equations.  Equations  (24)  then  become 


D  V  -  D  V 

S  S  «  « 


0-v.H-o 


(25a) 


(25b) 


and  In  similar  manner  the  normal  momentum  equation  and  the  energy 
equation  for  the  shock  structure  become 


DV 

s  s  3?  3c 


4  a_v 

3  a;2 


(25c) 


DV  |*  ‘i24  +  DVM|« 

SSd(  3.2  89 


i«  fw  3W. 
3C  1  3r 


-  0 


(25d) 


$fS8j '  ""  ‘V  w 


The  shock  structure  Is  taken  to  extend  from  the  outer  edge  of  the 
viscous  layer  ;  ■  0  to  the  freestream  ?  -*■  “.  At  the  shock  surface, 
the  shock  structure  flow  variables  and  their  normal  gradients  are  re¬ 
quired  to  match  the  corresponding  viscous  layer  values  at  n  *  n  as 

s 

given  by  the  solution  to  Eq.  (19). 

Following  Ref.  6,  we  now  illustrate  this  matching  procedure  at 
the  shock  surface  for  the  case  of  the  finite  tangential  shear  stress 
Ps(-|~)s.  In  terms  of  viscous  layer  variables,  this  stress  may  be 
written 


•  Mu)s  ‘  (5^s  •  <“ 


<“  FyV'a  -  <*!«>  Vt> 


where 


/[I  +  4  (cot  0  -  tan  0,)  sin  20,  ]  (tan  0  /cos  0,  )  for  j  ■  0 

I  L  8  b  D  8  D 


VJ)  * 


(1  +  (tan  0  / tan  0.)]  (3L/2x*) 

S  D 


for  j  *  1 


sin  0,  (1  +  tan  0  / tan  0,  )  (x*/L)  , 

W2U'  +  2 (tan  0  /tan  0,  ) ( 1  +  (tan  0  /2tan  0.)) 

s  b  s  b 


On  the  other  hand,  at  the  inner  edge  of  the  shock,  the  shear  is 


V 


. ***  . . . -  '''']ifliri'-'Ti||f|H 
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,  3u* 

(u  3/^s 


„2/3W. 

°-u-<Te)s 


(in_) 

v3y*;s 


(— ) 
v3u  'a 


lt2,3W. 

p-u-(a?s 


(30) 


with  the  subscript  “s"  here  denoting  conditions  at  r,  m  0.  Equat¬ 
ing  Eq.  (27)  to  Eq.  (30) 


Nf 

<TE>.  *  •*!«>  V”  A(lf>8  <31> 

This  matching  of  viscous  layer  solution  to  shock  structure 

solution  at  the  shock  surface  Is  shown  next  ior  the  case  of  the 

shock  surface  velocity  u*.  First,  recognizing  from  Fig.  2  that  at 

the  freestream  W  ■  cos  (6  +  0.  )  and  V  ■  -  sin (6  +  6,  ),  and  at 

00  s  b  ®  s  b 

the  shock  surface  W  ■  (u*/U  )  (u  /u*)  -  A  cos  0  ,  the  tangential 

S  S  o>  s  s  s 

momentum  equation  (25b)  is  Integrated  once  with  the  result 

aw 

(xr)  ■  [cos  (0  +  0  )  -  A  cos  0  ]  sin  (0  +  0.  )  (32) 

d;  s  s  b  s  s  b 


Substituting  Eq.  (32)  into  Eq.  (31)  and  solving  for  the  reduced 
velocity  A  at  the  shock  surface 


A 


u* 

s 

U 

09 


cos  0  2Q  (j)  Q  (j)  Nf 

/ - § - + - i - 2 -  / — mu  \  1 

cos  (0  +  0.  )  sin  [2(0  +  0,  )]  '  2f  * s 1 

S  D  s  b 


(33) 


For  the  case  of  flat  plate  flow  j  -  0  and  0.  *  0,  Eq.  (33)  reduces 

D 

to  the  result  of  Ref.  6. 

The  shock  angle  in  terms  of  the  density  and  pressure  behind 
tv<!  shock  is  obtained  next  by  integrating  the  normal  momentum 
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equation  (25c)  once  and  applying  freestream  boundary  conditions.  Neglect- 

3V 

ing  the  gradient  term  (— )  as  small  compared  to  V  [see  Ref.  5] 

sin  (0  +  0.  )  V  -  P  -  sin2  (0„  +  0.  )  -  (34) 

s  b  S  S  S  l)  ..t 

yM 

00 

Eliminating  V  by  means  of  Eq.  (25a),  the  shock  angle  is  then  given 
s 

by 


sin2  (0  +  0.  )  -  D  (P  -  l/yM2)/(D  -  1)  (35) 

s  b  s  s  •  s 

The  reduced  enthalpy  behind  the  shock  is  sought  next.  Integrat¬ 
ing  Eq.  (25b)  twice,  substituting  the  resulting  velocity  profile 
W(0  into  the  energy  equation  (25d),  and  integrating  twice,  we 
obtain  the  result 


s 


■  4  {[1  -  A  cos  0  /cos  (0  +  0.  )J2  cos2 

i  s  s  n 

+  sin2  (0  +  0.  )  +  2/(y  -  1)  M2} 
s  b  00 


<". +  V 


(36) 


where  the  gradient  term  (— )  has  been  neglected  compared  to  4> 

0  C,  S  S 

3W 

[see  Ref.  5],  and  where  the  term  (— )  has  been  eliminated  by  use 

p  £  8 

of  Eq.  (32). 

In  order  to  obtain  the  density  ratio  across  the  shock,  another 
2 

expression  for  sin  (0g  +  0^)  is  first  found  by  substituting  <!>g 


'--V 


to 


RMNMMM 
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from  Eq.  (36)  and  V  from  Eq.  (25a)  into  the  reduced  equation  of  state 

S 

P  -  D  (24»  -  V2)  (37) 

s  Zy  s  s  s 

Equating  the  resulting  expression  to  Eq.  (35) 

P8  ((Y  +  1)/(Y  -  D]  p8  +  i/ym2 

s  p«  Ps  +  (i  -  a  cos  eg/cos  (0g  +  6b)]2  +  (t/ym2)(y  +  1)/(y  -  l) 

(38) 

We  note  here  that  in  order  to  avoid  the  solution  of  a  cubic  equation, 

2  2 

a  term  D  [1  -  A  cos  0  /cos  (0  +  0,  )]  ,  which  appears  after  equating 

S  S  S  u 

2 

the  expressions  for  sin  (0  +  0.),  has  been  neglected  in  comparison 

S  D 

2 

with  (D  -  1).  This  is  justified  on  the  basis  that  since  the  re- 
s 

duced  shear  (Nf  )g  is  finite  but  always  small  compared  to  unity, 
then  by  Eq.  (33)  A  cos  0  -  W  is  always  close  to  W  ■  cos  (0  +  0,  ); 

S  S  ®  S  D 

and  the  larger  D  ,  the  closer  W  is  to  W  .  This  is  borne  out  by  the 

S  S  00 

numerical  results. 

The  reduced  pressure  P  [n  (x)]  at  the  shock  surface  has,  as 

s  s 

stated  before  in  conjunction  with  Eq.  (10),  been  assumed  equal  to 
the  reduced  pressure  P^Cx)  on  the  body  at  the  same  value  of  x. 

From  Eq.  (10)  written  in  the  form 

ye-  I(201/2  VW  I(\>  <39> 


and  from  the  definition  of  y  [Eq.  (2a)]  written  in  the  form 


-21- 


y8 


■I 


y  -  I  t(rK  +  y*  cos  0.  )/L]-*  dy*  -  { (r  +  r  )  x*  tan  0  /2Ly*]:y*  (40) 


w g  have  upon  elimination  of  y  and  the  use  of  Eq.  (14) 

s 


P„  ■  Pb  ■  -72  *  «1«>  «3<»  V1  <!> 
Pfl. 


(41) 


where 


(cot  6 


s 


Q3(j) 


2(x*/L)  cos  0. 


l3 C 2  +  (tan  0  /tan  0.  )]  (tan  0  /tan  0.) 

S  D  S  D 


for  j  ■  0 


for  j  *  1 


(42) 


and  where 


Tb  I(V 

T0  f(l1s> 


(43) 


Therefore,  if  from  Eq.  (19)  a  streamf unction  f(n)  and  its  derivatives 

were  to  be  known  at  n  ■  n  ,  the  relations  (33),  (35),  (38)  and  (41) 

s 

would  form  a  closed  system  of  four  equations  in  the  four  unknowns 

A,  0  ,  D  ,  and  P  . 
s  s  s 

IV.  Method  of  Calculation 

The  model  as  formulated  imposes  the  following  restrictions  on 
the  input  gas  and  body  properties:  M  sufficiently  large  such  that 

on 

2 

c  T.  %  U  /2,  T, /T.  «  1,  y  "  constant,  Pr  ■  1,  and  0,  <<  1. 

P  U  ®  DU  D 

The  solution  for  the  viscous  layer  streamfunction  f (n)  and  for 
the  various  flow  properties  of  interest  proceeds  by  first 
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transforming  Eq.  (19)  to  the  following  set  of  first  order  differential 
equations  as  in  Ref.  6: 

f  -  J  J  -  Z1/2(l  +  C.I)"1  g  g  -  -  Z1/2(l  +  C.I)'1  fg  (AA) 

n  n  z  n  2 

where 


Z  -  |l  -  J|  (1  +  C2J)  =  I 
The  no-slip  boundary  conditions 


(45) 


are  then  Imposed.  Now,  for  values  of  the  wall  shear  g(0)  •  (f  ), 

nn  b 

greater  than  the  value  which  gives  an  inviscid  flow  region  at  the 

outer  edge  of  a  boundary  layer  (n  -+•«),  a  velocity  profile  with 

finite  shear  g(n  )  ■  (Nf  )  at  the  shock  surface  can  be  obtained 
s  nn  s 

for  each  station  along  the  body  within  the  merged  layer.  At  each 

such  value  of  g(0),  the  no-slip  boundary  conditions  are  fully 

specified.  Beginning  then  at  n  ■  0  with  a  value  g(0)  •  [g(0)]+ 

n 

s 

only  slightly  larger  than  the  value  (g(0)]  which  gives  the 

n  -*•  °» 

9 

boundary  layer  solution  (n  -*■»),  Eqs.  (AA)  -  (A5)  are  simultaneously 

s 

integrated  out  from  the  body  surface  to  the  shock  surface  for  fixed 
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local  values  of  C^  and  C^.  The  Initial  guesses  for  these  parameters 

are  obtained  from  the  boundary  layer  limit  n  -►“with  0  guessed  to 

s  s 

be  8  *  0  where  0  ■  0  (y,  M  ,  0,  ,  j)  is  the  inviscid 

SS,  S,  S,  W  D 

inv  inv  inv 

shock  layer  value  for  wedges  and  cones  given  in  Ref.  12.  As  shown 
in  the  Appendix 


C 


1 


0 

S 


COS  (0 


[- 


cos  0 


_  +  e. ) 

inv  _ 2 


s 


inv 


0 

8  inv 


(47) 


0 

s 


0 

8  inv 


(48) 


The  numerical  integration  of  Eqs.  (44)  -  (45)  is  stopped  at  the  first 

point  n  reached  within  the  viscous  layer  where  J  ■  f  -  u*/u*  £  1. 

n  s 

The  coordinate  ng  is  the  finite  value  of  n  at  this  point,  and  the 

resulting  finite  values  of  f(n  ),  g(n  )  and  I(n  )  are  used  to  com- 

s  s  s 

pute  A  and  F  in  approximate  form  from  [see  Eq.  (33)] 


cos  0 


s 


inv 


2Qx(j)  Q2(j) 


(fi— )  ] 


-1 


cos  (0g  +  0fa)  sin  2(0g  +  0fa)  v2f'n"n 

* inv  inv  D  s 


(49) 


and  the  relation  given  by  Eq.  (43).  These  results  are  then  used  to 
provide  a  better  value  of  0g  from  Eqs.  (35),  (38)  and  (41).  Better 
values  for  the  local  parameters  and  are  next  obtained  from 
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the  relations  [see  Appendix] 


C 


1 


T0  COS  es  2VJ>  92<J)  >  -2 

t,  (cos  (e  +  e.  f  +  sin  [2(8  +  e  )]  *2f*n-n  1  <50) 

b  so  s  b  s 


[2  +  (tan  0  /tan  0,  )]  (tan  0  /tan  0.  ) 
_ s _ b _ s _ b_ 

2  "  l(n  )  f(n  ) 

s  s 


(51) 


and  beginning  again  at  n  *  0  with  g(0)  *  [ g (0) ]  ,  Eqs.  (44)  - 

s 

(45)  are  re-lntegrated  from  the  body  to  the  shock  surface  n  ■  ns* 
When  after  such  successive  Iterations  the  results  at  the  shock 
surface  from  one  cycle  agree  sufficiently  with  those  values  from 
the  previous  cycle,  the  local  no-sllp  viscous  layer  solution  Is 
available.  The  effects  of  small  local  velocity  slip  and  tem¬ 
perature  jump  are  then  computed  as  local  perturbations  from  the  no- 
sllp  results.  The  method  for  obtaining  the  local  slip  solution  Is 
entirely  analogous  to  that  given  in  Ref.  6  except  that  Eqs.  (47)  - 
(48)  of  that  reference  are  replaced  by 


u*,  sin  0.  . 

.JJ:  -  (f  )  r - r^-r  f  (1 

u*  un  o  l(3L/x*)J  lu 

s 


-  -^i)  Q 

3  ;  w3 


(j)f1  I 


jnr 


nsf, 


(Y  -  1)  T 


Vl/2 

T  * 


•  - ]n„  l(f  f)nnL.n 

2fnsl  Fnst  n  ns  nsl  nn  1)  0 


(52) 
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cos  6 


8nsi  ^nsi  .  1 

+  1 5 


cos  6 


s 


nsf,  .,-1 


si  cos  (6i  +  0  )  f  lA  .  cos  (0  +  0,)J 

s  .  b  nsf.  s  „  b  s 

nsH  nsf. 


to»,;  «» 


and  the  rarefaction  parameter  [Eq.  (35)  of  Ref.  6]  is  replaced  by 


M  C 
00 


1/2 


QX(J) 


Rex*!»  (1  ■  1A>  Vj)(Y  '  1)  F 


x*  sin  0.  , 


54) 


The  final  results  at  g(0)  ■  [ g (0) ]  m  constitute  the  local 

s 

solution  at  the  downstream  end  of  the  merged  layer.  The  final 
values  there  for  and  C ^  become  the  initial  guesses  for  these 
parameters  at  the  adjacent  station  a  small  distance  upstream.  The 
viscous  layer  equations  are  integrated  here  from  n  ■  0  with 
g(0)  -  [g(0)]+  +  Ag(0) ,  where  Ag(0)  is  a  small  increment  in 

n  w 

s 

the  surface  shear  g(0).  Through  this  stepwise  upstream  marching 

method,  with  successive  iteration  cycles  at  each  step,  the  velocity 

and  state  variables  are  computed  at  several  stations  x*  within  the 

merged  layer  regime.  The  computation  is  stopped  at  the  upstream  end 

when  the  condition  M  a  <<  1  or  the  condition  (f  ).  <<  1  is  no  longer 

°®  n  b 

satisfied.  If  the  computation  is  carried  to  stations  x*  farther 

upstream,  the  present  model  breaks  down  as  the  calculated  values  of 

p  /p  become  less  than  unity.  This  occurs  at  values  of  the  rarefac- 
s  ® 

1/2 

tion  parameter  M  (C/Re  .  )  of  the  order  of  unity. 

"  X"  , 00 
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VI.  Results  and  Discussion 

Computations  have  been  carried  out  in  the  manner  indicated  in  the 

previous  section.  The  cases  considered  cover  the  range  20  {  Ma  J  25, 

0.06  $  T./T  $  0.15  and  0°  $  6.  $  20°. 
d  0  b 

In  Figs.  3-5  the  density  at  the  shock  surface  for  sharp  wedges 

is  compared  with  measured  values  obtained  by  McCroskey,  Bogdonoff  and 
14 

Genchl  .  Considering  the  scatter  in  the  data,  the  agreement  with 
theory  would  appear  to  be  good.  The  theoretical  shock  density  ratio 
for  a  10°  cone  is  shown  compared  with  the  data  in  Fig.  6.  In  contrast 
with  a  wedge  having  the  same  flow  conditions  and  comparable  angle  6^, 
the  shock  at  the  same  distance  x*  from  the  leading  edge  is  weaker  on 
the  cone.  This  is  due  to  the  circumferential  spreading  of  the  flow 
around  the  cone,  bringing  the  shock  wave  closer  to  the  body  surface 
and  causing  the  merged  layer  to  extend  farther  downstream  than  on 
the  wedge. 

Fig.  7  presents  the  wedge  surface  pressures  given  by  Vidal  and 

9 

Uartz  in  comparison  with  the  theory  for  several  values  of  0,  .  The 

D 

strong  interaction  and  free  molecule  limits  given  by  Hayes  and 
13 

Probstein  are  also  shown.  Both  the  present  theory  and  the  data 
lie  below  the  result  predicted  by  strong  interaction  theory.  This  is 
also  seen  to  be  true  for  sharp  cones  as  shown  in  Fig.  8  where  the 


cone  surface  pressure  by  the  present  theory  and  by  the  strong  inter- 

19 

action  theory  of  Mirels  and  Elllnwood  are  compared  to  the  experi¬ 
mental  data  in  the  form  presented  by  Waldron^  and  by  Hof land  and 


Glick 


18. 


(A-2) 
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Using  Eqs .  (28)  and  (41)  -  (43)  to  eliminate  P  ,  we  obtain  then 

s 


C2  *  [2  +  (tan  0g/tan  0b)  ]  (tan  0g/tan  6b)/I(ns)  f(ng)  (A-3) 


from  Eq.  (20b). 

In  the  boundary  layer  limit  n  -*•  *  and  g(n  )  0,  we  find  that 

s  s 

since  f  'v  1  when  n  ^  n 
n  s 


[f(na)]|  -[  **.], 

1  r,  '■*  ®  i  1  n 

a  U  s 


(A-4) 


H(n5)] 


I  ’  ‘I 

1  r  -  ®  ' 


(1  -  f  )  (1  +  Cjf  )  dn] , 

n  'n  ♦  00 

0  S 


finite  (A-5) 


so  that  taking  0-0  ,  Eqs.  (A-l)  and  (A-3)  yield 

S  3inv 


s 


cos  (0  +  e, ) 

1  „  s.  b  „ 

0  , _ inv _ .2 

T  1  cos  6  1 

b  s. 

inv 


(A-6) 


0  "*0 
S  S  - 

inv 


s 

0  =0 
s  s . 

inv 


(A-7) 


Here  0  (y,  M  ,  0,  ,  j)  is  the  lnviscid  shock  layer  value  for  0 

S  .  «  D  8 

inv 

obtained  from  Ref.  12. 
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Fig.  1  Leading  edge  flowfield 
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Density  ratio  across  the  shock  structure  for  a  wedge 
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Density  ratio  across  the  shock  structure  for  a  cone 
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Fig.  11  Heat  transfer  rate  on  cones 


